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^ I A simple non-Hermitean random matrix (RM) model is used to study the Glasgow 
^ ' method of finite-density lattice QCD. The zeros of the RM partition function are evaluated 
through an averaging procedure, involving the zeros of the random 'propagator matrix' 
^ . in the complex chemical-potential plane. 

The nature of the uncertainty affecting the results is similar to that produced by round- 
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ing errors in computing the known analytic result. This similarity is exploited to give 
quantitative estimates on the relationship between the size of the matrix and the num- 
CN ' ber of configurations needed to achieve a given precision. For the quenched ensemble 
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considered here, the relationship is exponential. 



1. INTRODUCTION 



. Random matrix theory has had a significant impact on the theoretical and lattice study 
of QCD m]. This is especially true for lattice QCD at finite chemical potential 0. RM 
Qj \ models are solvable, easy to simulate, and reproduce some of the difficulties of lattice 
■ calculations. 

Lattice simulations at finite chemical potential n are plagued by the so-called sign 
^ . problem. The hermiticity of the /i-dependent piece of the Dirac operator p = If)^=Q + /i7o 
^ ! is opposite to that of the /i = part. This is due to the Euclidean formulation. As a 
result, the fermion determinant is in general a complex number and cannot be used as a 
probability measure. Therefore, generating an unquenched ensemble of configurations at 
nonzero chemical potential is not possible. 

The quenched approximation has been shown to fail in the context of the random matrix 
model employed here. Quenched simulations probe not QCD but the Nf — ^ limit 
of a theory where the fermion determinant is replaced by its absolute value [^. One 
may simulate the unquenched theory at nonzero /z, using an ensemble generated under 
different conditions, and including the fermion determinant in the observable. An elegant 
way to do this is provided by the Glasgow method [Q, with recent promising results 0. 

The unquenched QCD partition function at arbitrary yU can be seen as the expectation 
value of the fermion determinant, averaged over an ensemble weighted only with the gauge 
part of the action. One may also use an unquenched ensemble at zero chemical potential, 
but then the observable will also include the inverse fermion determinant at = 0: 

Z = (det(f(/i)))^„„^, = (det(^(/.))/det(|)(0)))^„„^,,^_^^^o (1) 
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The efficiency of the averaging process is determined by the overlap between the region 
of phase space probed by one's ensemble and that of the 'true', yU 7^ ensemble. The 
unquenched ensemble at = should do better in this respect than the quenched one. [] 

The idea of the Glasgow method is to expand the fermion determinant, which is an 
NxN matrix on the lattice, in powers of the fugacity ^ = {n enters in this combination). 
The expansion is finite and exact. It is obtained by rewriting the fermion determinant as 
det {p{fi)) = det (P + ^) where P is called the propagator matrix. The eigenvalues of 
P can be used to construct the expansion of det p. 

Since ^ is the same for all configurations, the averaged coefficients give the expansion of 
the partition function in powers of ^. Thus, with only one /i-independent ensemble, one 
has access to the full /x-dependence of the partition function. The location of the zeros in 
the complex /i plane maps out the phase structure. In particular, the zeros of Z{fj,) close 
to the real axis define the critical value(s) of the chemical potential 0. 

2. RANDOM MATRIX MODEL 

We consider a model defined by the partition function 

Z(m,fi) = JvCe-^^'^^^ det{D{m,fi)) ; D{m,fi) = 



m iC + jj 



(2) 



where C is an x complex matrix, and / PC stands for integration over each complex 
matrix element. This defines an ensemble of random matrices, which mimics the ensemble 
of gauge field configurations from QCD. The exponential factor ensures that the real and 
imaginary parts of each matrix element of C are independent Gaussian random numbers. 
D{m,ii) approximates the Dirac operator for quark mass m. It is anti-Hermitean for 
/i = 0, with no definite hermiticity otherwise. The chiral condensate and the number 
density are the logarithmic derivative of (0) with respect to m and n, respectively. 

2.1. The random propagator matrix and its eigenvalues 

Define the random analog, P{m) of the fermion propagator matrix as follows: 



P(m) = ^ det(D(m,/i)) = det(P(m) + yu). (3) 




Its main property is that its eigenvalues are the roots of the fermion determinant. Com- 
pute the baryon number density by substituting into (0) and taking the logarithmic 
derivative. The result is similar to the formula for the condensate using the Dirac matrix. 



Here {■ ■ ■)z stands for averaging with the partition function (^. The propagator matrix is 
not Hermitean in general, so its eigenvalues are scattered in the complex plane. For m = 0, 
they are uniformly distributed in the unit disc. Their distribution can be computed using 
the phase-quenched partition function, the same way it is done for the Dirac operator M. 



^It should be emphasized, though, that both approaches are in principle correct. Using the quenched 
ensemble is not equivalent to the quenched approximation, where the determinant is neglected altogether. 
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2.2. Glasgow averaging 

Consider an ensemble of N x N matrices C, generated with the Gaussian weight in 
For each C, we construct the corresponding fermion matrix and compute its eigenvalues. 
This way the fermion determinant in (0) can be written as a polynomial in fi, 

2N 

det(D(m,/i)) = det(P(m)+/i) = HiXk + fJ^) = I] Cfe/x^ (5) 

k k=l 

The average over the Gaussian ensemble performed on (|^) will commute with the finite 
sum. One obtains the partition function (0) as a polynomial in fi: 

2N 2N 

Z{m,^) = {det{D{m,fi))) Gauss = '^{ck)GaussfJ''' = Yl{Ak + fi). (6) 

k=l k=l 

Here, the stand for the zeros of the partition function. They may be interpreted as 
'average eigenvalues' of the propagator matrix. 

2.3. Exact zeros of the partition function 

The zeros of the partition function are accessible by direct calculation ^j^. The 
exact zeros for m = and = 96 are identical to those in the last panel of Fig.^ The 
zeros are located along the critical line Re{l + ^"^ + ln^) = 0. It can be derived from the 
formula, correct to lowest order in 1/N, 

Z(m u) = ^ —6-^^^+"'^ + - u2(^+i) (7) 

In the thermodynamic limit — > oo, the partition function (Q) splits naturally into two 
terms, given by simple analytic functions of the parameters. Each of them dominates ex- 
ponentially in some region of parameter space, corresponding to a thermodynamic phase. 
The derivatives of Z are discontinuous on the boundary separating these regions. In the 
thermodynamic limit, the zeros coalesce into a cut that follows this boundary. The critical 
value of fic = 0.527 (for m = 0) is determined by the point where the line crosses the real 
axis. 

3. SIMULATIONS 
3.1. Some results 

Results of Glasgow averaging for A^ = 16 and m = are shown in Fig.|I| . Complex NxN 
matrices C are generated with Gaussian weight. The eigenvalues of the corresponding 
propagator matrix are calculated. The polynomial whose roots are the eigenvalues is 
constructed for each 'configuration'. The coefficients are then averaged. The points 
shown in the figure are the roots of the averaged polynomial for a certain number of 
configurations. For comparison, the location of the exact zeros is indicated with crosses. 

The roots for a few configurations are very similar to quenched eigenvalues, since the 
eigenvalues of individual configurations are controlled by the quenched distribution. There 
is a clear indication that the result tends to evolve into the true zeros. There are two 
interesting points about how convergence is approached: (i) the cloud of zeros, even far 
from convergence, has a fairly regular shape; (ii) some zeros converge much faster than 
others. In our case, the zeros with the largest imaginary parts converge first. 
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Figure 1. Different stages of Glasgow averaging, for m = and = 16. 
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Figure 2. Effect of artificial noise on the exact zeros of the partition function. 

3.2. Artificial noise 

Since the coefficients are obtained through averaging, they are affected by a random 
uncertainty, similar to a rounding error. It is resonable to assume that the relative error 
is of the same order on each coefficient. One can study the effect of such an uncertainty 
by adding a fixed proportion of noise to the exact coefficients: Ck = Cfc(l + eRk) , where e 
is a fixed small number and Rk are uniform real random numbers between ±1. 

The effect of artificial noise is dramatic as seen in Fig.§ . A noise factor of the order 
of 10"^'' alters significantly the value of fic ■ These diagrams are strikingly similar to the 
last two in Fig.|l| . Again, the zeros closest to the real axis are the most sensitive. 

It is reasonable to expect that the effect of the errors is an exponential function of the 
matrix size N. For artificial noise, we can afford to check this. In Fig.^ the noise factor 
e has been 'tuned' to have approximately the same effect for each matrix size N, which 
is doubled from one frame to the next. The corresponding noise factor e is roughly the 
square of the previous one. Indeed, the artificial noise factor that has the same effect for 
different matrix sizes is proportional to a small number to the power iV : e^v ~ e~"^ . 
In other words, the precision in digits required to determine /ic with given accuracy, is 
proportional to the matrix size. 

3.3. Estimations and speculations 

Our assumptions need to be checked before making statements about simulations. Is the 
relative variance (j{ck)/ck of the coefficients always the same? From the existing results. 
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Figure 3. Artificial noise for different matrix sizes, tuned to lead to similar distortion 
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Figure 4. Simulations and artificial noise, so that the precision of Hc is the same. 



one can say that for given A^, the relative variance is fairly constant. For example, in a 
very long simulation for N = 16, it is very close to 4.00, except for the first few coefficients, 
for which it is smaller with a factor of 2.. 3. For larger values of A^, the relative variance 
tends in general to be constant within one set, and approximately given by log2 A^ . 

An important question refers to the correlations between coefficients. For instace, if all 
coefficients fluctuate proportionally, the effect on the roots of the polynomial is smaller 
than if the fluctuations are uncorrelated. At this point, there is no conclusive information 
on correlations. Assuming that they are negligible, we estimate the number of Glasgow 
configurations needed to get some fixed accuracy on /Zc , as a function of matrix size A^, 



log2 A^ Acfc 



-aN 



N. 



1/2 



conf 



Ck 



2 ^2aN 



(8) 



This estimate is hard to check, given the large numbers involved. Consider A^ = 16 and 
AT = 8. We have e,, = ej ^ {^/{N^ZY^') = {s/iN^H.Y/'f - A^i, 
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The comparison for artificial noise with eg = 0.07 and eie = 0.005, and simulations 
with 2000, respectively one million configurations is made in Fig.§ . Our estimations 
seem justified as all four panels have approximately the same error on /ic . For A^ = 256, 
which in terms of lattice size would correspond to a 4^ lattice with one degree of freedom 
per site, one gets €256 ~ 10""^° leading to a number of configurations around 10^° . 
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4. CONCLUSIONS 

In the context of a schematic random matrix model, it is found that the number of 
configurations needed to achieve a fixed precision using the Glasgow method, grows ex- 
ponentially with the matrix size A^. This is not surprising, since the ensemble of matrices 
used for the calculation is quenched. As a result, there is probably very little overlap 
between the region of 'phase space' probed by this ensemble and the one important in the 
unquenched ensemble at nonzero chemical potential. In terms of convergence, the proce- 
dure used here is equivalent to 'brute force' (using a quenched ensemble and including the 
complex determinant in the observable). The problem of achieving a large cancellation 
while averaging over a complex phase [§] is traded for the huge precision needed for the 
coefficients in order to obtain the zeros with reasonable accuracy. 

On the positive side, the Glasgow method is seen to converge here, albeit for small 
A^. The qualitative picture is slightly different from what is implied in The cloud 
of partition function zeros slowly narrows down to the true critical line, which is not 
distinguishable as long as the cloud is present. Some zeros converge faster than others. 
In general, the points in parameter space where the continuum limit partition function is 
larger are less sensitive to perturbations, leading to more accurate results. The limiting 
partition function is not zero along the critical line. The zeros along the 'higher' sections 
of the critical line converge first. 

An open question is to what extent an unquenched ensemble at zero chemical potential, 
like in lattice simulations, would improve the situation. The present model lacks certain 
symmetries that are present on the lattice. They may play a significant role, and can be 
included in a more sophisticated RM model. 
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